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UNEAR  AND  NONLINEAR  ELECTRON  CYCLOTRON  INTERACTION 
IN  OPEN  RESONATORS 

I.  INTRODUCTION 

One  of  the  most  recently  developed  devices  for  the  generation  of 
coherent  radiation  in  the  millimeter  and  centimeter  range  is  the  gyrotron 
(or,  electron  cyclotron  maser).  In  this  device,  radiation  power  is  pro¬ 
duced  by  the  interaction  of  the  TE  fields  of  the  cavity  with  an  electron  beam 
which  propagates  and  gyrates  along  an  external  static  magnetic  field.  This 
interaction  has  been  studied  extensively,  "  both  theoretically  and  experi¬ 
mentally,  with  impressive  results,  such  as  a  22%  efficiency  of  conversion  of 
beam  kinetic  power  to  radiation  at  a  2  mm  wavelength  with  22  kW  CW  output 
power.  Substantial  improvements  are  possible  if  axial  nonuniformities  are 

1 9  20 

introduced  in  the  form  of  gradients  either  in  the  oscillator  cross-section,  ’ 

21-23 

or  in  the  external  magnetic  field. 

The  frequency  of  the  gyrotron  is  determined  by  the  relativistic  gyro- 
frequency  of  the  electrons  in  the  external  magnetic  field.  For  efficient 
operation,  this  gyrofrequency ,  Doppler-shifted  in  the  electron  beam  rest 
frame,  must  be  in  close  synchronism  with  the  ei genfrequency  of  the  cavity. 
However,  since  the  size  of  the  cavity  cannot  be  made  arbitrarily  small,  a 
large  value  of  the  gyrofrequency  may  result  in  synchronism  with  more  than 

04  oo 

one  of  the  cavity  eigenfrequencies.  Although  preliminary  investigations 
have  not  been  conclusive,  it  is  expected  that  multi -mode  effects  may  prove  a 
problem  for  wavelengths  less  than  ~1  mm. 

To  avoid  such  potential  limitations,  it  is  proposed  that  the  gyrotron 
cavity  be  replaced  by  an  open  quasi-optical  resonator,  Fig.  1,  formed  by  two 
mirrors  of  approriate  shape.  Taking  the  z-axis  as  the  axis  of  the  radiation, 
the  fields  can  be  shown  to  be  derivable  from  the  usual  vector  potential 
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Fig.  1  —  Schematic  representation  of  the  quasi-opticai  resonator  and  electron  beam  configuration. 
The  diffracted  fields  are  contained  by  the  two  mirrors,  one  or  both  of  which  may  be  partially  trans¬ 
mitting.  The  individual  electrons  of  the  beam  execute,  to  zero  order,  a  helical  motion  in  the  exter¬ 
nal  magnetic  field  B0. 
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(1) 


A  =  2A(r,z)  sin  |kz  +  ct(r,z)][£x  cos(wt)  +  ee^  si n (cot )] 

+  ft7  fr  (A  cos(kz  +  a)]  [x  c°s(<*>t)  +  ey  sin(wt)]$z 

where  w  is  the  radiation  frequency,  k  =  w/c  is  the  dominant  term  of  the 
wavenumber,  a(r,z)  =  kr2/2R  +  tan“^(z/z0)  is  the  spatially  dependent  phase, 
A(r,z)  =  A0(r0/rs)  exP  (_'"2/rs)  is  the  position  dependent  amplitude  for 
the  fundamental  mode,  e  is  the  ell ipti city  of  the  radiation,  r  =  (x2  +  y2)2 
is  the  transverse  coordinate,  rs(z)  =  rQ(l  +  z2/z2)  is  the  sP°t  size  at 
z  with  rQ  its  minimum  value  at  z=0,  R(z)  =  z(l  +  z2/z2)  is  the  radius  of 
curvature  of  the  wave  front  at  z,  2Q  =  nr2/X  is  the  Rayleigh  length,  and 
X  =  2-rr/k  is  the  wavelength.  The  diffraction  angle  of  the  fundamental 
mode  is  0^  =  X/irr0.  The  diffraction  angle  of  the  higher  order  modes  is 
larger.  Therefore,  if  the  radius  of  the  mirrors  is  somewhat  larger  than 
rs,  most  of  the  energy  of  the  fundamental  mode  will  be  intercepted,  re¬ 
sulting  in  a  high  value  of  Q.  The  higher  order  modes  on  the  other  hand 
will  suffer  from  diffraction  losses,  their  Q's  will  be  small,  and  there¬ 
fore  these  modes  will  not  be  excited.  Thus  the  introduction  of  open 
resonators  provides  an  effective  natural  means  for  transverse  mode  selec¬ 
tion. 


This  paper  is  intended  as  a  preliminary  but  comprehensive  analysis 

of  the  electron  cyclotron  interaction  with  the  diffracted  fields  of 

Eq.  (1).  The  limit  we  have  considered  throughout  the  paper  is  that  of 

29 

an  electron  beam  propagating  along  the  radiation  axis.  We  have  assumed 
that  the  transverse  extent  of  the  beam  is  small  compared  to  the  radiation 
spot  size  r$,  therefore  the  action  of  the  field  on  the  beam  electrons  is 
that  of  a  plane  wave.  In  this  limit  we  have  developed  the  linear  theory 
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(Sec.  II)  and  the  nonlinear  theory  (Sec.  Ill)  of  the  interaction  in  a 
bounded  system  (oscillator).  Various  issues  pertaining  to  the  basic 
theory  are  discussed  in  Sec.  IV,  along  with  the  presentation  of  numerical 
examples. 

In  Section  II  the  linear  dispersion  relation  is  obtained  for  an 
oscillator.  For  the  sake  of  generality,  an  arbitrary  electron  beam  dis¬ 
tribution  is  assumed.  All  forward  and  backward  components  of  the  field 
perturbation  are  kept  and  the  value  of  the  refractive  index  n  =  kc/w  is 

left  unspecified  in  the  bulk  of  the  analysis  in  this  Section.  Specializing 

30 

to  the  case  of  a  monoenergetic  beam,  affected  only  by  the  forward  wave, 
yields  a  compact  result  and  reproduces  the  expressions  for  the  gyrotron 
in  the  cut-off  frame,  when  n  ■*  0.  For  luminous  waves,  n  -*•  1 ,  the  gain 
curve  has  two  interesting  features,  viz,  it  possesses  even  symmetry  about 
the  frequency  of  synchronism  and  can  lead  to  instability  even  at  exact 
synchronism,  provided  that  the  transverse  beam  velocity  satisfies  a  certain 
inequality. 

The  nonlinear  theory  is  developed  in  Section  III  for  a  constant 
amplitude  wave.  The  equations  of  motion  for  the  electrons  are  expressed 
using  the  axial  coordinate  z  as  the  independent  Lagrangian  variable. 

These  equations  are  interrelated  by  four  invariants.  Two  of  these  in¬ 
variants  relate  the  transverse  displacement  to  the  velocity,  while  the 
other  two  relate  the  velocity  components.  These  invariants  reproduce 
earlier  results  when  n  -*■  0  or  n  >  1.  A  transformation  to  slowly  varying 
variables  eliminates  an  additional  variable  and  results  in  a  single 
integrable  equation.  The  phase  space  plots  are  studied  and  distinction 
between  closed  and  open  orbits  is  made.  A  condition  for  complete  electron 


kinetic  energy  depletion  is  identified.  The  nonlinear  efficiency  is 
found  to  depend  on  only  four  combinations  of  the  field  and  electron  beam 
parameters  and  is  conveniently  presented  in  contour  plots. 

Finally,  certain  aspects  related  to  the  applicability  of  our  analysis 
are  discussed  in  Section  IV.  It  is  shown  that  the  symmetry  of  the  gain 
spectrum  can  be  used  to  relax  the  temperature  requirements.  The  interaction 
is  also  shown  to  have  inherent  abilities  for  axial  mode  selection.  In  our 
example  of  a  gyrotron  type  operation  we  obtained  a  13.7%  efficiency  at 
A  =  2  mm  with  a  61  kW  output  power.  The  potential  application  of  this 
interaction  as  a  free-electron  laser  is  discussed  and  a  scheme  for 
efficiency  enhancement  is  proposed.  By  utilizing  a  two-oscillator  con¬ 
figuration,  an  efficiency  n  =  35%  is  obtained,  while  estimates  predict 
efficiencies  can  reach  60%. 

The  cyclotron  interaction  of  electrons  with  luminous  waves  has  been 
the  subject  of  two  recent  investigations.31’32  In  the  first  paper  the 
dispersion  relation  was  obtained  and  analyzed  for  waves  with  arbitrary 
values  of  n  =  kc/w  in  the  limit  of  vanishing  axial  beam  velocity.  In  the 
second  paper  the  stability  was  investigated  in  terms  of  particle  energy  loss, 
obtained  from  an  expansion  of  the  equations  of  motion  in  the  radiation  field 
amplitude.  Both  these  analyses  assume  a  uniform  medium  of  infinite  extend  in 
z  and  use  the  time  as  the  independent  coordinate.  Accordingly,  their  results 
differ  from  ours,  since  in  the  present  work  the  axial  position  z  is  used  as 
the  independent  variable  and  the  presence  of  the  oscillator  boundaries  is 
important. 
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II.  LINEAR  THEORY 


The  fields  of  the  oscillator  are  given  in  terms  of  the  vector  po¬ 
tential  £  by  the  simplification  of  Eq.  (1), 


A  =  Aq  sin(kz)  (g.x  +  ie  ey)  +  c.c. 

(2) 

=  2Aq  sin(kz)  [|x  cos(wRt)  +  e  ey  sin(wRt) Je1 1, 

where  Aq  is  the  amplitude,  k  is  the  axial  wavenumber  (along  the  z  axis), 
iii  is  the  complex  frequency  consisting  of  the  real  frequency  and  the  growth 
rate,  w  =  wR  +  if,  e  characterizes  the  ellipticity  of  the  polarization, 
and  and  £  are  unit  transverse  vectors.  The  oscillator  fields  6JE  and 
6J3  are  given  by 


=  "  i  ft  =  1  c  Ao  sin(kz)  (£*  +  iy)  e'lwt  +  c.c., 

=  Z  *  A  =  -  ikA0  cos(kz)  (c|x+  i|y)  e"1wt  +  c.c. 


(3) 


In  addition,  a  static  uniform  axial  magnetic  field  is  present,  =  Bo£z> 

The  fields  given  by  Eq.  (3)  are  appropriate  for  those  in  an  oscillator 
formed  by  two  parallel  plates  of  infinite  extent  in  the  transverse  direc¬ 
tion.  They  are  also  applicable  to  a  general  case  of  an  oscillator  con¬ 
sisting  of  two  concave  boundaries,  provided  that  AQ  and  k  are  appropriate 
functions  of  x,  y,  and  z  and  the  small  axial  field  components  are  included. 
These  modifications  depend  on  the  waist  of  the  radiation  field,  r  ,  and 
the  Rayleigh  length  zQ  =  wr*/2c  and  are  negligible  if  the  beam  width  and 
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the  electron  L armor  radius  are  small  compared  to  r  ,  and  if  the  axial 
interaction  length  is  small  compared  to  zQ.  Accordingly,  under  such 
conditions  the  fields  in  Eq.  (3)  can  be  used  with  AQ  and  k  constant. 

The  behavior  of  the  field  is  governed  by  the  wave  equation. 


V2A  + 


'  9t2 


+  v 


k>*  = 


4tt 

c 


(4) 


where  v  is  a  real  quantity  representing  losses.  It  will  be  shown  later 
that  v  =  wr/Q,  where  Q  is  the  quality  factor  of  the  resonator.  Expressing 
the  current  density  as 

J  =  J  e”1u)t  +  c.c. ,  (5) 

it  is  found  that  Eq.  (4)  together  with  Eq.  (2)  can  be  put  in  the  form 


L 

2 


[( k2 c2  -  «J)  -  i(2T  +  v)wR]  Ao(ex  +  ieey) 


/ 


4ttc  J  sin(kz)dz. 


where  L  is  the  interaction  length. 


(6) 


Before  continuing,  it  is  necessary  to  evaluate  v  in  Eq.  (6).  This 
is  accomplished  by  integrating  Poynting’s  theorem  along  the  interaction 
region, 

L  L 

/ udz  +  Sz(z=L)  -  Sz(z=0)  =  -  f 6E-J  dz,  (7) 

O  0 
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where  u  is  the  radiation  energy  density  and  the  z  component  of  Poynting's 

vector  Sz  vanishes  at  z=0  and  has  the  value  Sz  =  (cjr/Q)  j udz  at  z=L,  by 

o 

the  definition  of  Q.  The  definition  of  r  implies  an  exponential  variation, 
exp(2rt),  for  the  energy,  thus 


(2P  +  Jp)  /udz  =  -  E  J  dz  , 


(8) 


and  from  a  comparison  with  Eq.  (6)  it  is  seen  that 


0)n 


(9) 


The  current  density  in  Eq.  (6)  will  be  evaluated  from  the  appropriate 
moment  of  the  distribution  function,  the  latter  obtained  from  the  solution 
of  the  Vlasov  equation.  Employing  the  expansion  f  =  f^  +  f^  of  the 
distribution  function,  f^  describes  the  zero-order  trajectories  in  the 
external  magnetic  field  and  is  therefore  a  function  of  the  corresponding 
constants  of  motion. 


p,2  =  p2  +  p2  , 

J.  x  y  ’ 

i  P„  iLm„ 

0  =  tan  (s^-)  — p —  z  , 
x  z 


P„  =  P,  . 


X_  =  X  -  w 

9  sz_m_ 

J  oo 


*g  =  y  +  Sjn 


o  o 


(10) 


which  are,  respectively,  the  transverse  momentum,  the  azimuthal  phase 
angle  relative  to  a  fixed  cyclotron  helix,  the  parallel  momentum,  and  the 
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x  and  y  locations  of  the  guiding  center,  defined  as  functions  of  the 

Cartesian  components  of  position  (x,y,z)  and  momentum  (p  ,  p  ,  p  ). 

x  y  z 

The  nonrelativistic  cyclotron  frequency  is  =  eB0/mQc.  Considering 
Eq.  (10)  as  a  transformation  from  Cartesian  to  constant-of-motion  phase 
space,  that  is,  defining  g^  =  f^  ang  g^  =  f^,  where  g^  and 
g^  are  functions  of  the  variables  defined  in  Eq.  (10)  as  well  as  z  and 
t,  the  Vlasov  equations  can  be  written  as 


(ft 


,(°) 


(11) 


Here  the  unperturbed  distribution  function  is  assumed  independent  of 
0,  Xg  and  yg.  The  relativistic  factor  y  is  defined  as  y2  =  1  + 

(Pi,  +  Pl)/(mQc)2. 

Using  Eqs.  (3)  and  (10)  and  introducing  the  expansion 


g^e- 


i  u>t 


+  c.c 


(12) 


in  Eq.  (11),  the  first  order  linear  homogeneous  differential  equation  for 
can  be  integrated,  yielding 


The  differential  operators,  D+  and  D  ,  acting  on  g^,  and  the  relative 


wavenumbers,  and  C2+»  are 


I  =  —  — —  +  — —  fp  I —  -  p  - —  \ 

+  c  3px  -  ymoc  yi  3p|(  9p 


=±k  + 


(nQ  -  yu))m0 


^2+  =  ±  k  + 


(fi0  +  Yw)mQ 


(14) 


The  constant  of  integration  in  Eq.  (13)  is  such  that  =  0  at  z=0, 
corresponding  to  an  unperturbed  entering  beam. 

The  perturbed  distribution  function,  obtained  in  Eq.  (13),  can  be  used 
to  evaluate  the  Fourier  component  of  the  current  density. 


J 


dp  . 


Performing  the  integration  over  e  in  Eq.  (15)  gives 


(15) 


(16) 

The  beam  density  n  has  been  extracted  for  convenience  from  the  distribution 
o 

function,  by  defining  nQgo  =  2irg^,  where  gQ  is  normalized  to  ydp1Jdp(fpAgo=i . 
Finally,  multiplying  j[  by  sin(kz)  and  integrating  over  the  length  of  inter¬ 
action,  yields  the  driving  term  required  in  Eq.  (6), 
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/ 


4tic  J  sin(kz)  dz 


"  *  T  eSVoc/dpj.dp/i 


(17) 

where  the  boundary  condition  sin(kL)  =  0  has  been  used  to  eliminate 
terms  similar  to  exp  [i(o1+  -  £-|_)L]. 

Equating  Eqs.  (6)  and  (17)  yields  a  consistency  condition  between 
r  and  e.  This  is  the  linear  dispersion  relation,  which  includes  the 
positive  and  negative  electron  cyclotron  interaction  with  both  the  for¬ 
ward  and  backward  traveling  waves.  However,  in  general  only  one  of  these 

interactions  is  resonant,  which  we  take  throughout  the  remainder  of  this 
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paper  to  be  with  the  forward  wave.  Defining 


f,  =  51+  -  k  + 


{%  -  Yw)mo 


P  =  o  i 


one  obtains  from  Eq.  (17), 


(18) 


ll 


L 


f  '!»,V  /dpidp«^  <£ 


%)  n 


+  e) 


(19) 


Comparing  this  equation  with  Eq.  (6),  it  is  seen  immediately  that  consis¬ 
tency  requires  that  e=l,  that  is,  only  a  circular  polarization  is  compatible 
with  the  interaction.  The  dispersion  relation  now  becomes 


U>r 


(u)£  -  k2c2)  +  iwR(2r  +  g~  ) 


*  16  P.  S  *  i2  • 


(20) 


where  cj2  =  4ue2n  /rn  and  u  is  a  complex  quantity,  defined  in  Eq.  (18)  and 
DO  0  0 

depending  on  u>  =  wR  +  ir. 


In  a  number  of  interesting  special  cases  only  the  real  part  uR  of 
p  is  of  significance.  These  cases  are:  (i)  the  (linear)  steady  state, 

T  =  0,  (ii)  the  conditions  just  above  threshold,  T  0,  and  (iii)  when  the 
beam  is  sufficiently  weak,  so  that  the  imaginary  part  of  u  is  negligible. 
In  such  cases,  separate  equations  can  be  obtained  for  the  frequency  and 
growth  rate. 


“R  -  (kC>!  ■  ?  moLC“bO  ^ 

“R<2r  *r|!s  '”oLc-u/dpj.dpu  ^ 


2yR  -  sin2uR 

(2uc)2 


sin2u 


"R 


(21) 


Finally,  for  the  special  case  of  a  monoenergetic  beam,  gQ  =  6(p^  -  p^Q) 
6(p a  -  Pn^/P^o’  integration  over  ^  and  p(|  gives 


12 


m 


4 


where  the  refractive  index,  n  =  kc/o>R,  has  been  introduced,  the  quantity 

uR  is  obtained  from  Eq.  (18)  in  terms  of  the  beam  quantities  yQ  and  p 

and  the  dimensionless  velocities  B(/0  =  P/JO/Y0moc  and  3X0  =  f>iphQ^Qc 

were  used.  The  operator  D  and  the  functions  R  and  Rn  are 

y  n 


To  complete  the  formalism  we  note  that  for  the  fields  of  Eq.  (3)  the 

total  steady  state  field  energy  inside  the  interaction  region,  L,  is 
-L 

given  by  J  udz  =  (ajRAQ/2c)2  L  (1  +  nz)/-n,  hence  the  efficiency  of  con- 
o 

verting  kinetic  energy  to  radiation  energy,  per  unit  time,  is 


We  have  thus  completed  the  linear  analysis  of  the  interaction. 

The  results  are  represented  by  Eqs.  (22)  and  (24)  for  a  monoenergetic 
beam  and  weak  (or  vanishing)  growth-rate,  as  well  as  the  general  case 
given  by  Eq.  (20).  Let  us  consider  now  the  significance  of  these  equations. 
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in  the  limiting  cases  n  ->  0  and  n  ■+  1 .  The  first  of  these  cases, 
n  -*  0,  corresponds  to  the  gyrotron  mechanism.  Referring  to  Eg.  (24), 
the  efficiency  of  the  oscillator  is  proportional  to  the  quantity 


Vq 


+  1  CL  /sinpR\2 

"R  )  2  C  %  \  ’^R  / 


(25) 


The  first  term  in  the  above  expression  is  always  negative,  and  its 
stabilizing  effects  are  reduced  if  the  ratio  B^0/8?0  is  large.  If  this 
ratio  is  large  enough  for  this  term  to  be  omitted,  a  positive  energy 
conversion  can  be  obtained  if  the  derivative  dR^/dv^is  positive.  This 
occurs  in  a  number  of  intervals  in  the  domain  of  pR  ,  the  most  prominent 
which  is  -  n  <  |jR<  0,  and  the  maximum  of  dRg/dp^=  0.54  is  attained  at 
uR=  -  1.3.  Such  negative  values  for  pRcorrespond  to  >  ^0/y0>  hence 
the  wave  frequency  must  exceed  the  relativistic  electron  gyrofrequency. 


Let  us  now  turn  our  attention  to  the  interaction  with  luminous  wave 
n  -+  1 .  In  this  case,  the  efficiency  of  the  resonator,  referring  to 
Eq.  (24),  is  proportional  to 


°mRQ 


Piio(1~g)io) 

®[o 


sinupV  d  /sinu^2 

Pr  /  d,R  ^  llR  j 


(26) 


This  function  is  plotted  in  Fig.  (2)  for  B  Q  =  0.1  and  (a)  B  Q  =  0.7, 

(b)  Bj0  =  0.425,  and  (c)  Bxo  =  0.35.  A  number  of  interesting  features 
can  be  observed  in  Eq.  (26).  First  it  is  noted  that  D  Rg  (i.e.,  n^) 
is  an  even  function  of  pR,  therefore,  contrary  to  the  gyrotron  mechanism 
the  interaction  with  luminous  waves  does  not  depend  on  the  sign  of  the 
frequency  mismatch,  AwQ  =  ojr  -  UQ/yQ  -  kv  =  -  2uRv(0/L.  Instability 
occurs  if  either  of  two  conditions  are  satisfied.  The  first  is 


NORMALIZED  LENGTH,  Mr 

Fig.  2  —  The  function  ),  which  is  proportional  to  the  linear  efficiency,  versus  *iR ,  which 

is  proportional  to  the  resonator  length,  for  n=l,  fJjj0  =  0.1  and  (a)  0io  =  0.70,  (b)  0|o  =  0.425,  and 
(c)/Jio  =  0.35. 


*lo  >  26  o  <'  -  eo>  • 


(27) 


and  yR  =  0.  If  Eq.  (27)  is  satisfied,  the  first  term  in  Eq.  (26)  is 
positive,  while  the  second  term  vanishes.  Thus,  contrary  to  the  gyrotron 
mechanism,  a  positive  efficiency  can  occur  for  Aw  =  0,  provided  that 
Eq.  (27)  is  satisfied.  If  Eq.  (27)  is  well  satisfied,  the  first  term 
dominates  up  to  the  point  |yR|  =  1.166,  where  tanyR  =  2  yR  and  therefore 
-  yR  dRq/dyR  =  Rg.  Additional  regions  of  instability  occur  when  the  second 
term  is  positive,  provided  that  Eq.  (27)  is  not  too  strongly  violated. 

These  regions  cover  the  intervals  jtt  <  ) y_ j  <  a.  -  1/a.,  where  a.  =  (2j+l) 

J  J  vJ 

ff/2,  j=l,2,3...,  and  the  upper  limit  is  a  good  approximation  to  the  solu¬ 
tion  of  tanyR  =  yR.  In  these  regions,  yR  dRg/dyR  reaches  the  maximum 

values  1/(1  +  b-)  at  |yD|  =  b.  -  l/(«/2  b.),  where  b.  =  (4j  +  1)tt/4,  and 
J  K  J  J  J 

j  is  a  positive  integer. 
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III.  NONLINEAR  THEORY 


The  linear  analysis  established  that  the  cyclotron  interaction 
with  plane  waves  results  in  instability,  even  for  luminous  waves,  n -*1 . 

In  this  section  the  nonlinear  equations  of  motion  are  obtained,  both  on 
the  fast  and  slow  spatial  scales.  These  equations  are  investigated  and 
constants  of  the  motion  are  obtained.  The  application  of  these  invari¬ 
ants  permits  the  analytical  study  of  the  conditions  for  energy  depletion, 
the  trajectories  in  phase  space  and  the  exact  scaling  relations  for  the 
efficiency.  A  single  first-order  differential  equation  is  found  to  be 
sufficient  to  describe  the  evolution  of  the  particles. 

The  nonlinear  equations  of  motion  are  studied  under  the  action  of 
a  uniform  axial  magnetic  field,  JBQ  =  BQ£2,  specified  by  its  nonrelati vistic 
electron  gyrofrequency ,  f2Q  =  eBQ/mc,  and  a  steady  state  wave,  whose  vector 
potential  is 

A  =  Aq  [fx  sin(kz  -  u>t  +  <J>Q)  +  e  |y  cos(kz  -  wt  +  4>0)J  .  (28) 

Only  the  forward  wave  is  included  in  Eq.  (28).  Adding  to  Eq.  (28)  a 
similar  expression,  with  k  -*■  -k  and  <j>  -*•  4>0  +  T*»  recovers  Eq.  (2).  However, 
the  backward  wave  has  been  omitted,  since  resonance  with  the  forward 
component  is  assumed. 

Although  the  natural  independent  variable  for  particle  motion  is 
time,  it  is  convenient  for  oscillator  configurations  to  adopt  a  Lagrangian 
representation,  expressing  the  equations  of  motion  in  terms  of  the  axial 
distance  z.  Since  dt  =  dz/vz,  the  phase  of  the  wave  as  a  function  of  z 
is  given  by 
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<Mz)  =  k 


-  syry) d2'  +  *0 


(29) 


where  Bz  =  vz/c.  The  equations  of  motion  become 


Aq  cosi^ 


dijj  % 

ciz  ■  r  bz 


d_ 

dz 


(yV  ■ 


-£ 


A  sintjk  4^-  +  — 


dz 


6, 


d_ 

dz 


(YSZ)  .  -k  A0 


cosiJj  -  e  -g-  sintli 

p,  i 


d_ 

dz 


Y 


cosiJj  -e 


(30) 


d_ 

dz 


x 


d_ 

dz 


.1 


5 


where  AQ  =  eAo/mc2,  Bx  =  vx/c,  By  =  vy/c  and  y2  =  (1  -  g2  -  e2  -  B*)_1 . 
Note  that  the  first  four  of  Eq.  (30)  are  interrelated  by  the  definition 
of  y.  In  addition,  the  equations  for  x  and  y  can  be  integrated,  in  con¬ 
junction  with  the  equations  for  yg  and  yBY,  respectively,  to  yield  the 
constants  of  the  motion 


lx  =  YBX  -  A0  sin*  +  ^y  , 


fl 

ly  =  YBy  "  eAo  C0Sli'  '  ~  x  • 

These  integrals  can  be  used  to  assess  the  assumption  of  transverse  uni¬ 
formity  implied  by  Eq.  (28). 
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An  additional  invariant  can  be  obtained  from  Eq.  (30)  by  noting  the 

30 

equations  for  yBz  and  y.  This  invariant  is 


i -I  *  y(n  -  ez)  , 


(31) 


and  its  invariance  is  associated  with  the  plane  wave  nature  of  the  fields. 

When  n=0,  this  invariant  describes  the  conservation  of  axial  momentum  in 

12 

the  cut-off  frame  of  the  gyrotron.  If  n  >  1,  I-j  is  proportional  to  the 

33  34 

conserved  particle  energy  in  the  wave  frame.  *  The  significance  of 
this  invariant  is  that  it  indicates  the  possibility  of  complete  energy 
depletion,  y=l ,  if  the  initial  particle  properties  are  such  that  I-j  =  n. 
Denoting  these  initial  properties  with  a  subscript  zero,  the  necessary 
condition  for  complete  energy  depletion  is 

-  n 

Yo  n-ezQ  .  (32) 


Introducing  the  definition  of  yq,  Eq.  (32)  can  be  written  as 


(' 


-  6 


zo 


n2+l  \ 

2n  )  ’ 


(33) 


where  B^Q  =  B2Q  +  B2Q.  An  interesting  observation  is  that  when  n=l, 

Eq.  (33)  reproduces  the  condition  for  the  vanishing  of  the  first  term  in 
the  linear  dispersion  relation  (see  Eq.  (26)).  It  should  be  noted  that 
Eq.  (32)  or  (33)  represents  only  a  necessary  condition  for  complete  energy 
depletion.  The  sufficient  condition  depends  on  whether  the  state  of 
complete  energy  depletion  is  accessible  by  the  particles.  We  will  return 
to  this  point  later. 
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The  introduction  of  the  invariants  reduces  the  system  of  Eq.  1 30 ) 

formally  to  two  coupled  integro-differential  equations  for  >2  and 

7 

Following  the  usual  procedure we  define 


=  ?,  s  i  n  ( t 1  +  x) 
x  e 

B  =  C0S(+  +  ,\) 


(34) 


where  B,  ,  v  and  x  are  functions  of  z  varying  slowly  over  a  cyclotron 
period.  Substituting  Eq.  (34)  into  Eq.  (30)  results  in  trigonometric 
functions  with  arguments  \  and  l.  >.  The  2.  +  y  terms  phase-mix  in 
a  few  wavelengths,  hence  the  slow  spatial  scale  equations  are 
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It  can  be  seen  that  when  n=0,  c=0  and  Sz  ->  0,  the  combination  1^  +  I '  +  1 

+  {Qq/u)z  gives  the  corresponding  gyrotron  type  invariant  in  the  cut-off 
12 

frame.  Also,  when  n  >  1,  then  transforms  in  the  wave  frame  to  the 
corresponding  second  invariant.  "^*34  pur^herrnore?  using  the  fast  length 
scale  Eq.  (30),  the  derivative  of  1^  is 


dL 

_ L 

dz 


(1-0  kA. 


,  o 

1 - —  (1  -  — ) 

nB,  '  -jjy; 


yB  sin(2y  +  v) 


,  1  i+f 

+  11  -  Hr}  ~  Vin2* 


(37) 


Hence  I2  is  invariant  also  on  the  fast  length  scale,  provided  e=l. 

The  determination  of  the  two  invariants,  I-j  and  I^,  allows  us  to 
study  the  trajectories  in  phase  space.  The  appropriate  coordinates  for 
the  two-dimensional  representation  are  x  and  y.  Using  the  definition  of 
y  and  the  invariance  of  1-j  and  Io  to  eliminate  Bz  and  in  terms  of  y  and 
X  (and  the  initial  conditions,  denoted  by  subscript  zero),  we  find  that 

C0SX  =  {]'n2)  (y~Yo)2  +  2(Aa>o/ai)  Yo  {y-yo]  +  (1+c)  Vo  B10  COSxo 
(1+e)  Ao  [(1-n2)  (y-y0)2  +  2(l-n320)  Yp  (y-yQ)  +  (y^)2]’2 

(38) 

where  Awq  =  w  -  UQ/yQ  -  kvzQ  . 

We  are  primarily  interested  in  the  case  of  luminous  waves,  n*l. 

In  this  case,  Eq.  (38)  takes  on  a  rather  simple  form.  The  resultant  equa¬ 
tion  is 

qCw  -  q2cosy  =  const  ,  (39) 
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where  Cw  is  a  constant  proportional  to  the  ratio  of  the  initial  fre¬ 
quency  mismatch  and  the  wave  amplitude. 


Yog-io 

>-ezo 


“  -  kvzo  -  no/To 
co(l+e)A0 


(40) 


and  q  is  a  linear  function  of  y. 


(y-Y0)  +  T. 


(41) 


It  can  be  seen  that  q=l  corresponds  to  the  initial  condition,  Y=Y0> 
hence  the  constant  in  Eq.  (39)  is  equal  to  Cw  -  cosxQ.  In  addition, 
depletion  of  the  perpendicular  or  of  the  parallel  energy,  correspond  to 
q=0  and  q=Cd  respectively,  where 


1  - 


28  (1-8  ) 
zov  zo ' 

62 

10 


(42) 


is  a  constant  characterizing  the  anisotropy  of  the  distribution.  When 
Cd=03  complete  depletion  of  the  energy  is  possible,  see  Eq.  (33). 


For  a  positive  frequency  mismatch,  that  is,  when  Cw  >  0,  the  phase 
space  trajectories,  obtained  from  Eq.  (39),  are  shown  schematically  in 
Fig.  (3).  The  arrows  indicate  the  direction  of  particle  motion.  Both 
open  and  closed  trajectories  are  possible  in  principle.  The  separatrix 
between  these  two  families  of  curves  is  shown  by  the  broken  bell -shaped 
curve,  described  by  the  equation  qC2  -  cos2x  =  0  and  defined  only  in  the 

interval  |x|  £  v/2.  The  largest  value  of  q  on  the  separatrix  is  at  x=0, 

_2 

q=Cw  .  The  vortex,  or  0-point,  of  the  closed  trajectories  is  at  x=0, 

Q  *  (2CW)'2. 
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RELATIVE  PHASE  ANGLE,  X 

Fig.  3  —  Trajectories  in  phase  space  coordinates  q,  The  dashed  curve  is  the  separatrix.  The  dashed 
straight  lines  give  the  initial  values,  q=l,  for  (a)  Cw>  1,  (b)  0.5  <  Cw<  1,  and  (c >  0  <  Cw  <  0.5. 


For  the  curves  in  Fig.  (3)  to  represent  actual  electron  trajectories, 

they  must  intersect  the  initial  condition  line,  q  =1.  The  location  of 

this  line  can  be  obtained  in  relation  to  the  value  of  C  .  If  C  >1, 

w  w 

then  the  peak  of  the  separatrix  lies  below  the  line  qQ=l  (line  (a)),  and 
all  particles  are  untrapped.  Trapped  particles  exist  if  Cw  <  1,  but 
always  more  than  half  of  the  initial  phase  angles  x0  correspond  to  un¬ 
trapped  particles.  Lines  (b)  and  (c)  in  Fig.  (3)  describe  the  initial 
conditions,  qQ=l ,  if  0.5  <  Cw  <  1  and  0  <  Cw  <  0.5,  that  is,  above  and 
below  the  vortex,  respectively. 

When  C  <  0,  the  phase  space  trajectories  are  obtained  from  those 
in  Fig.  (3)  by  a  reflection  about  the  line  x=71/2-  The  trajectories  for 
C  >  0  and  C„  <  0  are  shown  in  Figs.  4(a)  and  (b).  Because  of  this 
symmetry,  any  particle,  trapped  or  untrapped,  in  Fig.  4(a)  has  a 
corresponding  particle  in  Fig.  4(b),  which  executes  identical  changes  in 
q  (or  y).  This  recovers  and  extends  to  the  nonlinear  regime  the  linear 
result  that  the  evolution  is  independent  of  the  sign  of  the  frequency 
mismatch  (or  C  ).  The  phase  plot  in  Fig.  4(c)  corresponds  to  exact 
synchronism,  C  =0.  In  this  case  the  separatrix  reduces  to  the  lines 
X=ir/2  and  x=-t/2  and  all  particles  are  semi-trapped.  In  this  case  the 
relative  phase  x  is  bounded,  while  the  energy,  or  q,  is  bounded  only  from 
below.  Whether  the  initial  evolution  in  this  case  corresponds  to  net 
particle  energy  loss  or  gain  depends  on  whether  >  0  or  Cd  <  0,  as  has 

been  shown  in  linear  theory.  However,  it  is  clear  from  Fig.  4(c)  that 
eventually  all  particles,  regardless  of  their  initial  phase,  will  gain 
energy  and  will  be  tightly  bunched  in  x  (but  not  in  energy),  about 
X=-it/2.  Accordingly,  this  configuration  can  be  the  basis  for  a  particle 
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NORMALIZED  ENERGY,  q 


Cw>0 
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—n  0  n  —n  0  n  — n  0  n 

RELATIVE  PHASE  ANGLE,  X 

Fig.  4  —  Comparison  of  phase  space  trajectories  when  (a)  Cw  >  0,  (h)  Cw  <  0,  and  (c)  Cw  =  0. 
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accelerator.  The  details  of  such  an  application  are  currently  under 
study  and  will  be  discussed  in  a  separate  publication. 


The  determination  of  the  two  invariants,  I-j  and  Ig,  which  connect 
the  three  independent  momentum  components,  makes  possible  the  description 
of  the  electron  evolution  in  terms  of  only  one  quantity,  say  q  (i.e., 
the  energy).  When  the  index  of  refraction  is  unity,  i.e.,  n=l ,  the 
corresponding  differential  equation  is 


4a  = 


sgn(sinx) 


[q-(qcw  - 


cosxQ): 


"K 


q-c. 


(43) 


where  sgn  is  the  sign  function  and  the  reduced  length  c,  is  defined  as 


=  2Q+el.A  (LM  k ( z-z  )  . 

VX0  \  PjLO  / 


It  can  be  seen  that  Eq.  (43)  is  immediately  integrable.  The  inversion 
of  the  resulting  expression  is  straightforward  only  when  Cw=0,  Cd  <  0. 
This  case  is,  however,  uninteresting,  since  it  corresponds  to  energy  gain 
by  the  electrons. 

In  spite  of  the  difficulties  associated  with  the  inversion  of  the 
integral  of  Eq.  (43),  useful  information  can  still  be  obtained.  It  is 
noted  that  the  solution  of  Eq.  (43)  has  the  form 

q  =  q^max;  xo’  Cw’  Cd^  ’ 

where  cmax  is  the  value  of  c  corresponding  to  the  length  L  of  the  inter¬ 
action  region  in  Eq.  (44).  In  our  Lagrangian  representation,  the 

efficiency  is  given  by  n  =  - c  y-v0  )  /  (y0-1),  hence 
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(46) 


n^max’  Cw’  Cd’  V  =  Cn  » 

where  the  angular  brackets  denote  an  ensemble  average  over  x0  and  the 

coefficient  C  is  defined  by 

n 


c  = 


'o 


to 


ZO' 


(47) 


The  efficiency,  therefore,  does  not  depend  on  the  seven  parameters  k,A,e, 
L,  Bzo,  6  Q,  nQ  arbitrarily  but  only  on  cmax,  Cw,  Cd  and  Cn>  Further 
simplifications  arise  by  the  fact  that  is  just  a  multiplier  of  the 
efficiency  expression.  Large  values  of  are  possible  for  a  wide  range 
of  beam  properties,  as  can  be  seen  in  Fig.  (5).  In  this  figure  contour 
plots,  =  const,  are  presented.  In  addition.  Fig.  (6)  presents  contour 
plots  of  Cd  =  const. 


Before  proceeding  with  the  presentation  of  the  results  for  the  non¬ 
linear  efficiency,  it  is  worthwhile  to  rewrite  the  linear  efficiency, 

Eq.  (24),  in  terms  of  the  variables  defined  above.  We  use  the  definitions 

of  Cj,  C.,,  C  and  note  from  the  definition  of  c  that  it  is  related  to  u 
a’  vr  n  0 

by  =  -  2(1-Cj)u„»  where  un  is  defined  in  Eq.  (18)  in  terms  of  the 

W  iMa  X  U  O  U 

initial  values.  Then  from  Eq.  (24)  the  expression  for  the  linear  effi¬ 
ciency  is  obtained, 


n 


9, 


Cn 

T^cTI 


c2 

w 


<cd 


(48) 


This  expression  has  all  the  features  expected  from  Eq.  (46),  since  it 
involves  a  dependence  on  t  „  (through  u  ),  as  well  as  C,  Cj  and  C  . 

nia  X  O  w  u  n 
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0.5  1.0 

NORMALIZED  AXIAL  VELOCITY,  P*0 


Fig.  6  —  Contour  plots  for  the  anisotropy  constant  Cd  vs.  0^,  and  0|o 
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Let  us  proceed  now  to  solve  for  the  nonlinear  efficiency.  Two 
numerical  approaches  are  available,  the  first  approach  based  on  a 
numerical  integration  of  tne  complete  set  of  Eq.  (30)  to  obtain  the 
Lagrangian  particle  trajectories.  The  second  approach  requires  the 
integration  of  Eq.  (43)  only.  Both  approaches  have  been  used.  In  all 
cases,  the  agreement  has  been  excellent.  Equally  excellent  has  been  the 
agreement  of  the  results  with  linear  theory,  when  a  sufficiently  small 
amplitude  or  short  interaction  length  has  been  employed  in  the  simulations. 

The  dependence  of  the  normalized  efficiency  n/C^  =  l-<q>on  the  re¬ 
duced  distance  u0  is  given  by  the  solid  curve  in  Fig.  (7)  for  a  distribu¬ 
tion  characterized  by  =  -  2.5  and  a  wave  characterized  by  Cw  =  1.4. 
Comparing  this  curve  with  linear  theory  has  shown  that  initially,  up  to 
pQ  sl.6,  the  evolution  is  essentially  linear.  Beyond  that  point,  the 
efficiency  undergoes  oscillations  in  nQ,  reminiscent  of  the  corresponding 
linear  behavior,  except  that  these  nonlinear  oscillations  do  not  diverge, 
while  the  linear  result  of  Eq.  (48)  predicted  divergent  oscillations. 

The  location  of  the  relative  maxima  of  the  efficiency  vs^.  pQ  is  however, 
in  strikingly  good  agreement  with  linear  theory.  The  first  maximum  is 
at  yQ  =  3.7,  and  the  corresponding  optimized  nonlinear  efficiency  is 
n/C  =  0.164,  approximately  one  half  of  the  linear  value  np / C  =  0.300. 

For  these  parameters,  subsequent  maxima  correspond  to  negative  values 
for  the  efficiency.  The  dashed  curve  corresponds  to  the  values  =  -2.5 

and  C  =  2.7.  This  case  will  be  discussed  later. 

w 

The  contour  plots  of  Fig.  (8)  show  the  normalized  nonlinear 
efficiency,  1  -  <q>,  at  the  first  maximum  (the  one  near  uQ  =  3.7),  for 
various  values  of  the  parameters  and  Cw.  As  can  be  seen  from  the 
Figure,  large  normalized  efficiencies,  e.g.,  1  -  <q>  ■  0.15,  are  possible 
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NORMALIZED  LENGTH,  Mo 

Fig.  7  —  Spatial  evolution  of  the  normalized  efficiency  l-<q>  =  rj/Cn  for  Cd  =  -2.5  and 
Cw  =  1.4  (solid  curve),  as  well  as  for  CH  =  -2.5  and  Cw  =  2.7  (dashed  curve). 


BEAM  DISTRIBUTION  CONSTANT,  Cd 

Fig.  8  —  Contour  plots  of  the  normalized  efficiency  l-<q>  vs.  the  parameters  Cd  and  Cw 
for  an  interaction  length  corresponding  to  the  first  efficiency  maximum  (p0  =*  3.7).  The 
dotted  curve  gives  the  value  of  Cw  for  maximum  efficiency  at  any  value  of  Cd. 


over  a  wide  range  of  beam  distributions  and  ratios  of  frequency  mismatch 
to  radiation  field  amplitude.  In  addition,  the  dotted  line  in  Fig.  (8) 
gives  the  value  of  Cw  required  for  optimal  operation  at  a  given  value  of 
C^.  For  such  values  of  Cw>  the  normalized  efficiency  as  a  function  of 
Cj  is  given  in  Fig.  (9).  The  overall  maximum,  1  -<q>=  0.165,  is 
obtained  for  =  -  2.2,  Cw  =  1.5. 

The  results  presented  so  far  apply  to  the  first  maximum  of  the  non¬ 
linear  efficiency,  at  a  point  near  yQ  =  3.7.  This  is  not  meant 
to  imply  that  the  first  efficiency  maximum  is  always  dominant.  For 
example,  referring  again  to  Fig.  (7),  it  can  be  seen  that  increasing  the 
value  of  Cw,  causes  the  amplitude  of  the  first  efficiency  maximum  to 
decrease,  and  that  of  subsequent  maxima  to  increase.  When  the  point 

Cw  =  2.7  is  reached  (dashed  curve),  for  which  1  - < q>  maximizes  near 
w 

yQ  =  6.8,  the  normalized  nonlinear  efficiency  is  equal  to  1  -  <q  >  =0.143. 
For  even  higher  values  of  Cw,  it  is  the  third,  fourth,  etc.  efficiency 
maximum  that  dominates.  The  same  behavior  is  observed  if  -  is 
increased,  keeping  C,.  constant.  Contour  plots  for  the  efficiency  at  the 
second  and  subsequent  maxima  are  very  similar  to  those  of  Fig.  (8), 
except  that  the  curves  are  shifted  to  slightly  higher  values  of  Cw  and 
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OPTIMAL  NORMALIZED 
EFFICIENCY,  1  -<  q  > 


BEAM  DISTRIBUTION  CONSTANT,  Cd 


Fig.  9  —  Optimal  normalized  efficiency  (i.e.,  along  the  dotted  curve  of  Fig.  (8)),  vs.  Cd. 
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IV.  DISCUSSION 


In  the  previous  two  Sections,  we  have  considered  the  linear  and 
nonlinear  electron  cyclotron  interaction  with  (primarily)  luminous 
waves.  From  that  analysis,  results  have  been  obtained  regarding  the 
stability  of  these  waves  and  the  nonlinear  steady-state  performance. 

In  this  Section,  we  conclude  with  various  aspects  concerning  the 
practical  feasibility  of  a  device  based  on  this  analysis.  We  show  that 
this  interaction  has  weak  temperature  requirements  (6Bzo/ezo  <<  ])  and 
that  excellent  mode  selection  can  be  achieved.  One  example  of  a  cyclo¬ 
tron  maser  is  given  and  the  use  of  this  interaction  as  a  free-electron 
laser  is  discussed.  Finally,  the  possibility  of  improving  the  nonlinear 
efficiency  (up  to  n  ^0.60)  is  discussed. 


Let  us  consider  first  the  effects  associated  with  a  thermal 


spread,  66ZQ  and  6yQ,  of  the  beam  electrons.  The  sharpness  of  the 
resonance  in  both  linear  and  nonlinear  theory  is  characterized  by  the 
parameter  uQ,  which  can  be  written  as 


p 


o 


(49) 


where  X  =  2-rr/k  is  the  radiation  wavelength.  From  the  dependence  of  the 
efficiency  on  yQ,  it  is  seen  that  |  6;iq |  «  -n/4  is  required.  With  regard 
to  the  parallel  velocity,  this  translates  to  the  requirement 


(50) 
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For  conventional  cyclotron  maser  operation,  w  >  f2/y0  +  kcBz0  and  this  require¬ 
ment  gives  63Z0/6Z0  «  A/4L,  prohibiting1^  operation  in  a  system  longer  than 
a  few  wavelengths.  However,  when  to  =  kc,  the  syrrmetry  about  yQ  =  0  allows  a 
choice  of  operating  point  with  go  <  PJyQ  +  kcBZQ.  In  particular,  if  the 
wave  is  suppressed  (e.g.  ,  by  a  polygonal  ray-path  in  a  multi-mirror  resonator) 
and  the  choice  is  made  that 


uo  =  11  (X  +  4}  »  (51 ) 

where  a  local  maximum  of  the  efficiency  is  located,  the  requirement  is 
simply 


(52) 


(53) 


which  may  be  more  restrictive. 

As  regards  the  other  significant  parameters  of  the  interaction,  the 
requirement  «  1  -  gives  a  condition  supplementary  to  Eq.  (53), 


(54) 


while  the  condition  on  the  parallel  velocity  is  less  restrictive  than  (52). 
Similarly,  the  condition  |6CW|  «  Cw  is  satisfied  if  the  above  conditions 
are  satisfied,  except  that  in  the  highly  relativistic  limit,  Bzq  -*•  1 ,  an 
additional  condition  is  imposed,  viz. 


(55) 
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(b)  Mode  Selectivity. 

Transverse  and  longitudinal  mode  selectivity  can  be  dealt 
separately.  We  will  not  deal  in  any  detail  with  the  former.  It  is 
expected  that  by  an  appropriate  design  of  the  mirrors,  only  the  lowest 
order  transverse  mode  will  be  intercepted,  while  the  higher  order  modes 
will  suffer  from  diffraction  losses. 


As  far  as  the  longitudinal  modes  are  concerned,  we  first  seek  to 
determine  the  number  of  modes  that  fill  the  positive  gain  curve  segment, 
at  whose  peak  the  desired  longitudinal  mode  is  located.  Since  the 
half-width  of  the  gain  curve  is  |6y|  &  it/4,  it  can  be  seen  from  Eq.  (49) 
that  the  number  6£  of  longitudinal  mode  pairs  (in  addition  to  the  mode 
at  which  operation  is  desired)  within  the  positive  gain  segment  is 


6£  < 


Bzo 

zn^y 


(56) 


Thus,  no  competing  axial  modes  are  present  if  f?zo  <  2/3.  On  the  other 
hand,  many  modes  will  be  excited  in  the  case  of  a  highly  relativistic 
electron  beam  with  8Z0  -*•  1. 

An  additional  point  to  be  considered  is  whether  the  operation  of 
the  device  at  a  certain  efficiency  maximum  is  accompanied  by  the  presence 
of  a  different  longitudinal  mode  which  might  correspond  to  a  value  of 
yQ  lying  within  a  different  positive  segment  of  the  gain  curve.  This 
question  is  particularly  important  when  there  are  many  wavelengths  within 
the  resonator.  In  that  case,  the  operation  is  seen  from  Eq.  (51)  to 
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require  a  large  value  of  uQ  and  the  broken  line  in  Fig.  (7)  indicates 
that  more  than  one  gain  segment  corresponds  to  positive  gain.  Such 
undesired  modes  can  be  suppressed  by  the  introduction  of  distributed 
feedback.  Inherent  means  for  mode  selection  are  also  available.  Thus, 
it  can  be  seen  from  Eq.  (49)  that  if  the  choice 


6zo  =  2j  +  1 


(57) 


is  made,  where  j  is  an  integer,  then  the  closest  potentially  unstable 
mode  falls  at  the  next  ( 2 j - 1 ) - th  maximum  of  the  gain  curve.  If  j  is  not 
too  small,  then  it  is  expected  that  this  maximum  may  correspond  to  nega¬ 
tive  gain,  thus  guaranteeing  the  absence  of  any  competing  modes. 


( c )  E xample  of  cycl otron  maser  operation. 


Before  proceeding,  it  is  necessary  to  present  a  numerical 
example,  in  order  to  demonstrate  how  our  results  can  be  applied  to  a 
particular  design.  We  consider  an  interaction  at  the  third  maximum 
of  the  gain  curve.  As  can  be  seen  from  Eq.  (51),  temperature  effects 
are  reduced  by  choosing  a  system  length  equal  to  3\.  The  normalized 


efficiency  at  the  third  maximum  is  n/C^  =  -.156,  for  =  -  2.9, 

Cw  =  +  3.0  and  pQ  =  +  10.1.  A  value  Ezo  =  0.143  (=1/7)  guarantees  the 
absence  of  any  competing  modes  (see  Eq.  (57))  and  Eq.  (53)  gives  the 


condition  6yq/y0  <<  1.2%.  From  the  definition  of  Cd,  Eq.  (42)  (or,  from 
Fig.  (6)),  it  can  be  found  that  B^0  =  0.25,  hence  yq  -  1-0433  correspond¬ 
ing  to  a  beam  energy  of  22.2  keV.  Then  from  the  definition  of  C  , 

Eq.  (47),  or,  from  Fig.  (5),  it  is  found  that  Cn  =  0.88.  Therefore, 
the  nonlinear  efficiency  is  n=  0.137.  For  the  above  values  and  the 


definition  of  Cw  (Eq.  (40)),  the  normalized  vector  potential  of  the 
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radiation  field,  assuming  linear  polarization,  has  the  value 

O 

Aq  =  1.55  x  10  .  The  external  magnetic  field  gives  a  gyrofrequency 

satisfying  /w  =  1.055  and  the  frequency  mismatch  is  AcoQ/u) 

=  (oj  -  fi  /y  -  kv  )/w  =  -  0.154.  Finally,  by  choosing  a  wavelength 
A  =  2  mm,  this  example  requires  an  external  magnetic  field  of  Bq  =  56.5 
kG  for  a  radiation  field  amplitude  of  EQ  =  833  sV/cm  (=  250  kV/cm) 
and  frequency  oj  =  0.942  x  1012  sec'^  (=  150  GHz).  For  an  incident  beam 
current  of  20  A,  the  radiated  power  is  61  kW. 


where  EQ  =  ojAq/c  is  the  amplitude  of  the  electric  component  of  the  radi 
tion  field.  The  operating  frequency  is  essentially  given  by  the  first 
term  in  Eq.  (58).  Defining  the  axial  relativistic  factor,  Yzo=0-B20)~ 
the  operating  frequency  is  accurately  given  by 


a)  =  (1  +  8  )y2  — 

v  zo' rzo  yq 


and  bears  the  same  form  as  the  free-electron  laser  frequency, 
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employing  a  magnetic  wiggler. 
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Since  tne  relativistic  gyrofrequency  depends  on  the  mass  factor 


Y  »  the  actual  quantity  that  multiplies  the  nonrelativistic  gyrofrequency 

to  give  the  operating  frequency  is  Ij"',  where  I ^  =  Y0(1~B20)  is  the  first 

invariant.  For  high  operating  frequencies,  it  is  required  that  »  1. 

This  requirement,  together  with  the  desire  that  be  not  too  much  less 

than  unity  (in  order  to  achieve  adequate  efficiencies)  requires  that 

-Cd  «  1.  Thus,  as  can  be  seen  from  Fig.  (9),  efficiencies  of  the  order 

of  10%  are  to  be  expected.  For  example,  w/Qo  =  10  requires  I-j  =  0.1. 

Taking  yQ  =  50,  gives  £Zq  =  0.998  and  B  =  0.06.  Then  Cd  =  -0.11,  giving 

o/Cn  =  0.102,  and  =  0.917,  and  the  result  is  an  efficiency  n  =  9.36%. 

(Of  course,  higher  efficiencies  can  be  obtained  by  adopting  efficiency 

enhancement  schemes,  as  discussed  in  IV(e)  below.)  In  this  example,  it 

is  required  that  C„  =  3.6  and  u  =  3.6.  Choosing  a  wavelength  of  A=0.2  rrni 
w  o 

(f  =  1.5  THz)  and  a  system  length  L  =  lm  =  5000  X,  it  is  found  that  an 
external  magnetic  field  of  strength  BQ  =  47.5  kGs  is  needed.  Furthermore, 
assuming  a  beam  current  of  10  A,  the  radiated  power  is  equal  to  Prad  =23.4 
MW. 


( e )  Efficiency  enhancement. 

The  nonlinear  analysis  of  Sec.  Ill  shows  that  the  performance 
of  an  oscillator  operating  with  luminous  waves  is  limited  to  efficiencies 
n  &  16%.  Here  we  will  outline  briefly  a  scheme  which  can  improve  the 
performance  of  the  interaction  dramatically.  We  introduce  a  two-stage 
oscillator  configuration,  in  which  the  first  stage  (the  prebuncher)  acts 
to  pre-bunch  the  particles  in  their  relative  angle  x»  without  significantly 
changing  particle  momentum.  The  second  stage  (the  converter)  converts 
the  kinetic  energy  of  the  bunch  to  radiation  energy.  A  100%  conversion 
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is  possible,  if  the  bunch  is  perfect  and  =  0. 

In  the  prebuncher  we  require  a  combination  of  a  "large”  frequency 
mismatch  and  "small"  radiation  field  amplitude,  such  that  |C  -j|  >>  1. 

In  this  stage  the  electrons  are  highly  untrapped  and  only  negligible 
energy  exchange  occurs.  To  first  order  x  is 


where  =  0  has  been  assumed.  For  x0  uniformly  distributed  in  the 
interval  (0,  2tt),  Eq.  (60)  describes  the  formation  of  a  bunch. 

The  extraction  of  particle  energy  occurs  in  the  converter.  From  the 
discussion  of  Sec.  Ill  (see  also  Fig.  (3)),  it  is  recalled  that  a  particle 
on  the  separatrix  eventually  loses  all  its  kinetic  energy,  if  Crf  =  0,  as 
is  assumed  here.  Since  the  bunch  has  a  spread  in  x  rather  than  in  q,  it 
is  necessary  to  arrange  the  converter  parameters  so  that  Cw2  =  1.  With 
an  appropriate  choice  of  Cw-| ,  the  bunch  will  be  centered  at  the  peak  of 
the  separatrix.  The  center  of  the  bunch  will  thus  be  able  to  convert  all 
its  energy  to  radiation  energy,  by  choosing  the  appropriate  length  for 
the  converter  section.  Examining  the  bunch  coherence  from  Eq.  (60)  it  can 
be  estimated  that  efficiencies  of  the  order  of  60%  are  attainable. 


As  an  example  of  a  two-stage  oscillator  we  consider  the  case 
Cw]  =  7.5,  Cd  =  0,  with  a  discontinuity  in  the  radiation  field  amplitude, 
cw2  =  i.o,  for  c  >  The  evolution  of  the  efficiency  for  this  case  is 
shown  in  Fig.  (10).  The  oscillatory  curve  gives  the  evolution  when  the 
first  oscillator  is  sufficiently  long,  c-|  00-  However,  if  the  location 

of  the  discontinuity  is  appropriately  chosen,  the  ensuing  evolution  is 


NORMALIZED  EFFICIENCY,  1  -<q> 


Fig.  10  —  Efficiency  enhancement  in  a  two-stage  oscillator  configuration.  The  oscillatory 
curve  gives  the  nonlinear  normalized  efficiency  vs.  normalized  distance  for  Cwl  =  7.5  and 
Cd  =  0  in  a  single  resonator.  The  introduction  of  a  second  resonator  with  Cw2  =  1.0  at 
?=?! ,  shown  by  the  small  circles,  improves  the  overall  efficiency.  The  cases  shown  corre¬ 
spond  to  (a)  ?!  =  0.90,  (b)  ?i  =  1.72,  and  (c)  ?x  =  2.50. 


42 


characterized  by  substantially  increased  efficiency.  Thus,  for  the 
choices  ^  =  0.90,  1.72,  and  2.50,  shown  in  the  figure  by  the  small 
circles,  the  evolution  of  the  efficiency  in  the  second  oscillator,  where 
Cw2  =  ^*0*  is  shown  by  the  curves  (a),  (b),  and  (c).  This  procedure 
yields  the  optimal  values  of  the  efficiency  n  =  0.15,  0.28,  and  0.35  at 
the  positions  ^  =  1.5,  2.3,  and  3.1,  respectively.  This  clearly  demon¬ 
strates  the  principle  of  efficiency  enhancement  in  a  two-stage  oscillator. 
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